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Abstract. We numerically study the classical and quantum dynamics of an atomic bright 
soliton in a highly-elongated one-dimensional harmonic trap with a Gaussian barrier In the 
regime of the recent experiment by Dyke et ai, the system realizes a coherent nonlinear 
soliton beam-splitter and interferometer whose accuracy we analyze. In the case of tighter 
radial trap confinement and enhanced quantum fluctuations, a split soliton can represent a 
spin-squeezed, or alternatively, a fragmented condensate with reduced phase-coherence that 
can be measured by interfering the split soliton by the barrier. We also find large quantum 
mechanical uncertainties in the soliton's position and momentum due to nonlinear interaction 
with the barrier. 
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1. Introduction 

Bright solitons are solitary waves robust to collisions with each other and with external 
potentials, and are commonly studied within many physical systems including nonlinear 
optics [II10 and Bose-Einstein condensates (BEC) JH ffl [H [6l|. In the BECs bright solitons 
form peaks in atomic density against a negligible background density and are found in 
condensates of attractively-interacting atoms, or repulsively-interacting atoms in an optical 
lattice (band gap solitons) 0. In two-species BECs a bright solitary wave can be stabilized in 
one component while trapped inside a co-propagating dark solitary wave (phase kink) in the 
other component, forming a dark-bright vector soliton [l8ll9l [T0l[m[T2l[T31 . 

Despite the considerable theoretical and experimental efforts, the main research emphasis 
on bright solitons has been in classical mean-field properties, e.g., within the Gross-Pitaevskii 
equation (GPE). Soliton dynamics in BEC experiments (31 H [5]| have been adequately 
described using the GPE, although quantum and thermal fluctuations are suggested as a 
seed for condensate collapse mechanisms [[T4l| and soliton-train formation [15J. Recently, 
however, it has been suggested that quantum fluctuations cause the fragmentation of a soliton 
train [fT6l [TTIl and may already be more important in the present experiments than previously 
thought. Quantum effects in bright soliton dynamics have been studied in a nonlinear optics 
context [fT8l[T9ll20ll2ni22ll23ll24| . and also in few atom strongly-interacting systems (far from 
the regime we consider in this paper), specifically macroscopic superposition (Schrodinger 
cat) states of solitons [|25ll26l . tunneling through a barrier in a Bose-Hubbard model ETll . and 
soliton entanglement [|28ll . Quantum effects have also been predicted in gap solitons If29l[30ll. 

Motivated by the recent bright soliton experiments [31] we study the dynamics, 
reflection, transmission, and splitting of bright solitons in a combined trap of a highly- 
elongated, prolate harmonic potential and a Gaussian optical barrier. The experimental 
conditions of Ref. [i3Tl correspond to large atom numbers in a one-dimensional trap which is 
not radially very tightly confined. Quantum fluctuations in the ground state of such a system 
are weak and we first analyze the classical mean-field dynamics in this limit. Displacing the 
condensate from the equilibrium position at the trap centre initiates dipolar oscillations of the 
BEC which forms a propagating soliton. We introduce a Gaussian barrier at the trap centre 
over which a soliton is either completely transmitted, reflected or split into a transmitted and 
reflected component. The system realizes a strongly nonlinear, phase-sensitive soliton beam- 
splitter and interferometer. It has been previously suggested that the robust character of bright 
solitons lends itself to precision measurements, e.g., of surface reflectivity [l32ll and for atom 
interferometry ||33]| . 

In order to analyze the atom number statistics of the soliton interferometer formed 
by the laser barrier acting as a beam-splitter, we extend the theoretical study beyond the 
classical mean-field approach by studying the soliton dynamics in the stochastic phase-space 
representation within the truncated Wigner approximation (TWA) 041 [3511361 [371 [381 [39ll40ll . 
Although the output of the highly nonlinear beam-splitter is very sensitive to the relative phase 
difference between the solitons suggesting a possibility for high precision interferometry, the 
numerical analysis of the atom statistics in the TWA simulations shows that the nonlinear 
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soliton-barrier interactions enhance fluctuations of the relative soliton atom numbers in 
the interferometer output above those given by the classical limit of binomial statistics. 
This enhancement of fluctuations is large enough to cancel any benefit from the improved 
sensitivity to the phase values and happens even in the limit of negligible phase decoherence. 

We also show that, in a radially more tightly-confined trap and with smaller atom 
numbers, quantum fluctuations can be significantly enhanced, resulting in dissipative soliton 
dynamics. From the TWA simulations we calculate ensemble averages that provide quantum 
mechanical expectation values and uncertainties of the soliton's position and momentum, as 
well as the relative atom number statistics and phase-coherence of a soliton pair. Moreover, 
individual stochastic realizations of the TWA dynamics represent possible outcomes of single 
experimental runs, revealing considerable variations within an ensemble, e.g., with one 
trajectory displaying transmission of a soliton through the barrier and another one reflection. 
We also find that, due to the interaction with the barrier, when the soliton's kinetic energy 
dominates the nonlinear interaction energy the soliton may split into a spin-squeezed soliton 
pair propagating to opposite directions. For enhanced fluctuations the initial condensate 
undergoes fragmentation into a mutually incoherent soliton pair. We show that the loss of 
phase-coherence in the condensate fragmentation can be measured by studying atom number 
fluctuations in the recombination process of the solitons by the barrier. Finally, we also 
propose an approach that could potentially exploit quantum phenomena in the spin- squeezing 
of the solitons and overcome the shot-noise limit of interferometers based on classical atomic 
samples. 

2. Soliton Experiments 

Experiments on atomic BECs with attractively interacting atoms have produced single or 
multiple bright solitons, using Li-7 JH [3l |3T]|, or Rb-85 [5] atoms trapped in a cigar-shaped 
geometry. In all of these experiments, a Feshbach resonance was used to tune the s-wave 
scattering length of the inter-particle interactions to the desired negative value, producing a 
controllable attractive interaction between the atoms. A stable condensate takes the form of 
a single soliton [3, 31 J. Otherwise, if the attractive interactions are too large to be stabilized 
by the quantum pressure of the trapped atoms, the condensate collapses, leaving behind a 
predictable number of solitons which may each oscillate within the trap and collide with 
each other many times during an experimental run [[5]|. In the experiments in Refs. [|4l[5l| the 
stability of multiple solitons was attributed to the relative-phase of n between neighbouring 
solitons that leads to repulsive interactions between the adjacent solitons and the characteristic 
collective oscillatory dynamics of the solitons [|T5l . The essential condition for the well- 
defined relative phase of n to remain throughout the dynamics is that the solitons were 
assumed to be sufficiently phase-coherent. The phase-coherence of solitons has been the 
subject of some previous theoretical investigations [,14J , and in particular it was suggested in 
Ref. [16J that the loss of phase-coherence in the soliton trains could cause rapid fragmentation 
of the condensate. 

The recent experiment OTTl . which forms the starting point of our theoretical study, is an 
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ideal system to study the soliton phase-coherence. In this experiment, a BEC comprising a 
single soliton of 2 x 10^ ^Li atoms was formed in the centre of an elongated ID harmonic trap 
that was generated by a tightly-focused laser beam. The BEC was made weakly attractive by 
tuning the scattering length to a ^ -0.5ao via a Feshbach resonance, where ao denotes the 
Bohr radius. The radial and axial angular frequencies were cOr = 2;rx300Hz and co = 2;rx3Hz, 
respectively, such that although the system was not strictly ID, it was stabilized against 
collapse, with the chemical potential /j. < ficOr, resulting in significantly reduced density 
fluctuations along the radial direction. The potential minimum along the axial direction was 
suddenly displaced such that the soliton was placed in a nonequilibrium position, inducing 
dipolar oscillations for the condensate in the trap. A thin barrier, generated by a near-resonant 
cylindrically focused laser beam was inserted at the centre of the trap. The propagating soliton 
interacted with the barrier, and depending on its kinetic energy and the size of the barrier, 
underwent a reflection, transmission, or was split into two. The split solitons propagated in 
opposite directions before re-colliding at the barrier where they were observed to output in 
either direction [demonstrated in the numerical simulations in Fig.[Ua)-(c)]. In Sees. [3] and [5] 
we investigate the splitting behavior and the phase-dependence of the solitons' re-collision at 
the barrier first using classical mean-field analysis and then going beyond the mean field by 
including atom statistics. 

The soliton-barrier system provides an interesting atom-optical realization of a coherent, 
nonlinear soliton beam-splitter and interferometer. The nonlinear classical dynamics in this 
system can also lead to the formation of a bound pair of solitons or a 'soliton molecule' as 
a result of the interaction of a soliton with the barrier [41 J. An alternative method to split 
a single soliton coherently into two by means of applying density modulation with resonant 
pulses between difi'erent internal states was proposed in Ref. [|42||. 

3. Classical Dynamics 

We first investigate the classical periodic dynamics and splitting of a soliton in a harmonic 
potential within the GPE. The system may be considered to be in a ID regime when the 
chemical potential fi and thermal energy ksT are small compared with the transverse trap 
energy hco,-. For simplicity, here we assume that the system is sufficiently tightly confined 
in the radial direction so that the radial degrees of freedom are confined to a Gaussian radial 
ground-state wavefunction, and the dynamics are described by the ID GPE: 



where g = Itiojra, tOr is the angular frequency of the radial trapping potential, Vext is the 
external axial potential including the axial harmonic potential V = mo/x^/2, the Gaussian 
barrier [Eq. ([T5l) 1 and, in some simulations, a phase imprinting potential [Eq. (|2T])1. 

In the absence of the external trapping potential, the ID GPE with attractive interactions 
(g < 0) supports soliton solutions with the centre-of-mass position for the soliton j denoted 
by qj and its centre-of-mass velocity vj = dgj/dt. Using the notation x = x/l,-, qj = qj/h, a = 




2m dx^ 



(1) 



Quantum dynamics of atomic bright solitons 5 




Figure 1. (a)-(c), (f)-(i) Atomic density in the ID harmonic trap as a function of time for 
a propagating soliton within the classical mean-field analysis in the absence of quantum 
fluctuations. The soliton is initially displaced from the trap minimum and hits a potential 
barrier at the trap centre, subsequently splitting into two solitons that propagate into opposite 
directions. The two solitons then collide at the barrier at the trap centre and, depending on 
their relative phase, emerge either to the negative x direction, positive x direction, or the both 
directions. In (a)-(c) the nonlinearity takes the value \ U\ - 3fiLo, the barrier width cr - 0.021/ 
and the barrier heights Vo for the splitting and solitons' re-collision are both close to lAAtio). 
In (f)-(i) \U\ - \5tLU), cr - 0.14Z, the barrier height Vo - Tifico for the splitting, and the final 
barrier height Vq - lATiu). The kinetic energy per atom in the soliton as it collides with the 
barrier is Te - llfiu). The right-hand soliton is multiplied by a complex phase: (a) = 0, (b) 
= 1.2 (c) = 2.6, (f) = 0, (g) = 2.0, (h) = 2.4 and (i) = 2.7. The diff'erent stages 
of the interferometer scheme are marked on (b). (d) Relative atom number of the outgoing 
solitons 5N scaled by twice the total atom number 2N as a function of the phase 4> applied by 
direct multiplication of the soliton wavefunction before the second collision with the barrier 
(red line with triangles) for nonlinearity \U\ - Jifito, and by ramping a constant potential 
across one of the solitons [see Eqs. ( l20b and ( l2Tl il for \U\ - itio) (green line with squares) 
and \U\ - IQtico (blue line with crosses), and for comparison the exact solution (cos(f>)/2 
[see Eq. (l24b l for U - (black line) with numerical points from the GPE simulation (black 
circles) for the corresponding parameters, (e) as (d), but for a thicker barrier with cr = 0.14/. 
The trajectories in (a)-(c) have phases chosen such that at the end of the dynamics, all the 
atoms populate the left/right outgoing soliton or are evenly distributed between both outgoing 
sohtons. The trajectory in (f) has the imprinted phase = 0, and the trajectories in (g)-(i) span 
the region in (e) with large gradient, and have imprinted phases as described above. 
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a/lr,T = t(x)r,Vj = dqj/dr, where = (JhlmajrY^^, the soliton solution may be written in the 
following form Il43ll4^ : 



\a\N'^ 

^j{x, t) = ^j-^sech [\a\Njix - qj)] (T^^-'^'V'^^'^ , (2) 



with 



0y(T) =A^JaV + 5/T) + 0f, (3) 

Sj{T) = Tv]l2. (4) 

Here Nj is the number of atoms in the jth soliton, and Yij^j = where is the total 
atom number in the system. The spatially independent part of each soliton's phase has been 
encapsulated in 0,(r) and 0^°^ denotes the initial value of the phase. For the case of a single 
initial soliton, the initial phase value 0^"^ represents the global phase of the wavefunction and 
is inconsequential to the dynamics. If several solitons are present, they may approach each 
other, collide, and then leave each other unscathed up to shifts in position and phase. In the 
harmonically trapped system we consider, we calculate the ground-state soliton profile by 
imaginary-time propagation of the GPE, and find it to be slightly distorted compared with 
the homogeneous solution. In the presence of a harmonic trap the dynamics is no longer 
integrable, but it is easy to show that in a weak harmonic potential, for the approximate 
solution of a displaced soliton, qj and vy oscillate harmonically, such that 

qj = qjQ cos {ojt) + Vjq sin {ojt) (5) 
Vj = vyo COS (cur) -I- qjQ sin {ojt) (6) 



where qjQ and Vjq are the initial position and velocity of the jth soliton, co = oj/iOr, and the 
term, 

Sjir) = 



2a) 2 



sin(2a>r)^^^^^^2a,r)-^, (7) 



2 2^2 

replaces Eq. @i in the phase dynamics of Eqs. ([3]) and (O [45]. From the soliton wavefunction 
solution given by Eq. Q we can calculate the expectation values for the first and second 
powers of soliton position and momentum as 

(Xj), = J x\ilfj{x)\'^€j^dx = qj, (8) 
(xjh = J x'l^jixfej^dx = + q], (9) 

r dij/j{x) _, 

{Pj)s = -ifi I lAyW dx = mvj, (10) 

{pjh = -n' J ^j{xr^!^e-:'dx = + m^vj, (11) 



where 



ej= I \ilij{,xfdx = Nj. (12) 
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Here, the notation {...)^ indicates we are averaging across one soliton wavefunction. These 
relations provide the corresponding position and momentum uncertainties of a single soliton 
in the classical GPE mean-field solutions 

(Axj), = Ux/), - {xj)^V" = , (13) 

^ ' l^aNj 

(Apj)s = ({pj\-{pj)l) =-^. (14) 

These uncertainties are determined by the width of the soliton's wavefunction in Eq. ©. 

We model the experiment in which a single soliton was split into two solitons of 
approximately equal size (in this paper, by soliton size we refer to the number of atoms in 
the soliton). We assume the atoms are initially in equilibrium in the centre of the trapping 
potential and are described by Eqs. Q and (l5l|7]). We numerically simulate the displacement 
of the soliton by a distance Xd (equivalent to the experimental trap displacement of -Xd), 
followed by the application of a Gaussian barrier of width cr and height Vq, 

n(x) = Voexp|-^j, (15) 

at the trap centre. The soliton starts propagating from the edge of the trap towards the 
barrier. We investigate the behavior of the colliding soliton by varying the ratio between 
the nonlinearity \U\ and the centre-of-mass kinetic energy per atom, Te, for the soliton as it 
hits the barrier, defined by 



U=?^. l=J±. (16) 
/ y mo) 

TE=^-mo?xl. (17) 

Here / denotes the axial harmonic oscillator length and Xd the initial centre-of-mass 
displacement of the soliton from the trap centre. The expression (flTl) represents the initial 
potential energy when the centre-of-mass velocity of the soliton is zero. This potential energy 
is assumed to be converter to kinetic energy at the trap minimum. 

We study the soliton dynamics using parameters close to those of the recent experiment 
on . For example, in the experiment the Gaussian laser barrier of a 5//m waist (1/e^ intensity 
radius) in the trap with a» = 2;rx3Hz corresponds to the width cr =^ 0.16/ in Eq. (fT5l) . A typical 
experimentally studied value for the ratio between the kinetic energy of the soliton and the 
nonlinearity was Te/\U\ - 5.5 (for \U\ - AShto and the initial trap displacement of 0.5mm, 
or Xd - 231). We vary the barrier width between the very narrow values of cr ^ 0.012/ 
or cr ^ 0.021/ and broad values with cr ^ 0.14/ or, in Sec. 15. 4[ with cr ^ 0.32/. For 
the narrow barrier case the interaction time of the soliton with the potential is minimized, 
representing a sharp ideal barrier, reminiscent of a (5-function potential. Here cr =^ 0.14/ is 
close to the experimental values of Ref. OTTl . For ^Li atoms in the trap of Ref. OTll with 
oj = In X 3Hz, the narrow barrier width of cr 0.021/ can be achieved with a diffraction- 
limited focusing of a laser with the wavelength 650nm or shorter. In a weaker trap the barrier 
may be experimentally realized with a less tightly-focused laser. In most simulations we vary 
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the nonlinearity and use the initial displacement - 5.91, such that the kinetic energy as the 
soliton hits the barrier Te - llUco. For the nonlinearities \ U\ - 3hco this yields approximately 
the same ratio TeI\U\ - 5.5 as in Ref. Bill . For the largest nonlinearity we examine 
{\U\ = ISfioj) we also consider smaller displacements = 3.21 such that Te - S.lhco. 

When the nonlinear energy dominates (we examine the case I i7 1 = 15?Z(X>andr£ ^ 5.2%oj) 
the soliton does not split but may become stuck to the barrier for a finite time, and/or be either 
completely transmitted or completely reflected gaining a njl phase-shift on reflection [i.e., 
the reflected soliton is obtained by replacing (p^^^ by 0^°' + n/l and vy by -Vj in Eq. ©J. The 
njl phase shift may first appear surprising, but it is analogous to the phase difference between 
a reflected and transmitted electromagnetic wave from a finite-width dielectric slab enclosed 
between two optically identical media [|46l . 

Lowering \U\ or increasing the initial displacement (as in the example case we studied 
with Te - llfioS) the soliton may split into two solitons (a transmitted part and a reflected part, 
the latter of which gains a njl phase-shift), as weaker attractive interactions are less inclined to 
hold the soliton together. The proportion of the reflected and transmitted component depends 
on the interactions, as well as the barrier parameters (see Fig. |2l). We define the relative atom 
number ^A'^ and relative phase 6(f) between the solitons by 



where A^2 and (p2 are the number of atoms and the phase [see Eqs. Q and ([3])] in the left soliton 
(x < 0), and A^i and (p\ are the number of atoms and the phase in the right soliton {x> 0). 

We display the splitting dynamics as a function of the interactions and the barrier height 
and width [Vq and cr in Eq. ([T5l) 1 in Fig.[2l In the specific parameter regime shown in Fig. [2a), 
for a given barrier height a larger proportion of soliton is transmitted for small \U\ than for 
large \ U\. Consequently, a larger barrier is required to split a soliton into two equal parts when 
the interactions are reduced. Figure [2tb) demonstrates the intuitive result that a barrier of a 
smaller width has to be higher to split a soliton into two equal parts. 

After a soliton is split into two solitons of equal size, which subsequently evolve in the 
harmonic trap, the solitons re-collide at the barrier with a relative phase of 6(f> = n/l. At 
this point, the solitons are half the size of the original soliton and will interact in the same 
manner but as if the nonlinearity \U\ were half the strength. To eliminate this effect, we 
reduce the barrier height between the splitting of the first soliton and the subsequent collision 
of the solitons at the barrier to a height that would split a soliton of size N/l into two equal 
parts. If the solitons undergo free evolution in the harmonic trap between the splitting and 
recombination, they usually recombine to form a single soliton [see Fig.[Ua)]. We also find, 
however, an exceptional case, where the nonlinearity is large, and the barrier is sufficiently 
wide for a prolonged nonlinear interaction with the soliton [e.g., the trajectory in Figinf) with 
parameters a = 0.14/ and \U\ = ISfico]. In this case, the solitons do not fully combine, but 
two solitons emerge from the collision at the barrier. Moreover, in the case of wide barriers a 
soliton can collide inelastically, such that even for a barrier chosen to split a soliton into two 
parts of equal atom number the system can be asymmetric. In this case the reflected soliton 



6N = N2-Nu 

6(f> = 02 - 01, 



(18) 
(19) 
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Figure 2. (a) Relative atom number between two solitons 6N in the classical mean-field 
dynamics. The two solitons are formed by splitting a single soliton into transmitted and 
reflected components by a potential barrier, defined in Eq. dTsl l. The relative atom number is 
normalized by twice the total atom number 2N and is displayed as a function of the potential 
barrier height Vo- The width of the barrier is o" = 0.14Z and the kinetic energy per atom of the 
soliton as it hits the barrier is Te - llTicj. The nonlinearity (occurring in Eq. ([1])) is |t/| - 3tia) 
(squares), 5fia> (triangles) (circles) and IS/zw (crosses), (b) as (a), but with \ U\ = 3fiLL>, 
and with a barrier of width cr ^ 0.021/ (squares) and cr ^ 0.012/ (circles). 



is slower than the transmitted soliton, and performs lower amplitude oscillations [see Fig. 

m-m 

We study the relationship between the relative phase Sep and relative atom number 6N in 
the outgoing solitons by applying a constant relative phase difference by directly multiplying 
the wavefunction of one of the solitons by the phase factor exp [icfi], such that its wavefunction 
i/zj becomes tffj exp [icp] . We find that generally, the solitons do not fully recombine, but 
there are two outgoing solitons of different sizes with 6N that depends on the imprinted 
phase. Hence, we can use the system as an interferometer in detecting weak forces that 
accumulate a phase difference between the solitons. Since the two split solitons are spatially 
separated, in principle they can experience a difference in some weak force, e.g., in gravity 
that generates the change in the relative phase between the solitons. The accumulated relative 
phase difference can be measured by observing the relative atom number between the reflected 
and transmitted parts after the collision with the barrier. In order to simulate such a phase 
imprinting procedure which may be performed in an experiment we consider applying an 
example potential of the form 



V = Vr{t)@(x), 



where 



Vr(t) = y„ exp [-fi^it - tmf] , 



(20) 



(21) 



between times t = and tf, where t^ = {tf + ^/)/2, and @{x) is the Heaviside step function, so 
the potential is felt by one soliton and not the other; classically, in the absence of any phase 
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fluctuations, this potential imprints the phase 



Jti 



(22) 



where Erf is the error function defined by: 



2 

Erf(z) = — I exp \-y^\ dy . 



(23) 



'n Jo 



We choose /3 = co and show in Fig. [TJd) that the phase imprinting by the example potential 
yields almost identical results to multiplying one soliton wavefunction by exp [i(p]. 

The dependence of the relative atom number on the relative phase between the two 
solitons is easiest to understand in the case of a non-interacting gas (i.e., U = 0). We assume 
that in the first collision the barrier splits a soliton of containing A'^ atoms into two equal 
parts and in the second collision the barrier is changed so that it splits a soliton containing 
N/2 atoms into two equal parts. If we assume that the phase (f> is imprinted on the initially 
reflected soliton, that soliton has then acquired a phase difference -6(p = + n/2 with respect 
to the initially transmitted soliton. After the second collision with the barrier, the soliton 
propagating to the negative x direction consists of a transmitted component of the initially 
reflected soliton and the reflected component of the initially transmitted soliton (which now 
acquires a n/l phase shift upon the reflection). Adding both contributions together results in 
the atom number N cos^(0/2) for the soliton propagating to the negative x direction, where N 
denotes the atom number of the initial soliton. The atom number of the soliton propagating to 
the positive x direction is sin^(0/2) and the relative atom number after the second collision 
with the barrier reads 



We find that in the presence of nonlinear interactions the dependence of the splitting on 
the relative phase diff"ers from the noninteracting result. Moreover, the splitting behaviour 
is different for different barrier widths. If the solitons are split and recombined by a 
narrow barrier (we simulate a barrier of width cr =^ 0.021/), the dependence of the relative 
atom number in the output solitons on the phase difference is qualitatively similar to the 
noninteracting case, but the sensitivity to the imprinted phase is enhanced. As shown in 
Fig. [Ud), when the nonlinearity \U\ is increased, the functional dependence deviates from 
the exact sinusoidal behaviour of the non-interacting case due to the increasingly nonlinear 
interaction with the barrier. As | U\ increases, the value of ^ = 0^ at which 6N/N is a maximum 
is reduced from (p = n at U = 0. Consequently, the gradient of 6N/N becomes steeper for 
< (/),„ and shallower for > 0„,. Further narrowing the barrier has only a small effect on 
this deviation: using a barrier of approximately half the width (reduced from cr =^ 0.021/ 
to cr - 0.012/) and height Vq - lOOfito (chosen to split the soliton into two equal parts) 
only reduces the maximum difference in 6N/N between the \U\ = 3tico and U = cases 
by 3%. For a wider barrier (we simulate a barrier of width cr = 0.14/ that is close to the 
experimental values of Ref. OTll ). at weak nonlinearities, e.g., \U\ = 3fico, the deviation 
from the noninteracting result is similar to that of a narrow barrier case. However, as the 




(24) 
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nonlinearity is increased, e.g, to\U\ = ISfico, the curve [5N/N](4>) becomes strongly distorted 
from a sinusoid. We also find that the first minimum in SN/N is shifted from = to - 2, 
as shown in Fig.[Ue). The significant asymmetry of the curve representing the relative atom 
number as a function of the imprinted phase is reminiscent of the nonlinear behaviour of 
bistable systems. This could indicate the existence of more than one solution of ^A'^ for a 
given value of (f>. 

The steep gradient of 6N/N as a function of the applied phase for the nonlinear 
solitons indicates that measuring the relative atom number between the solitons exiting the 
barrier could potentially provide a very accurate detection method for the applied phase and, 
hence, an accurate weak force detection. The classical analysis of the mean-field dynamics 
based on the GPE simulations, however, does not yield any information about the statistical 
fluctuations of the atom numbers, resulting from the nonlinear interactions of the soli ton with 
the barrier. Additionally, quantum fluctuations may increase phase fluctuations between the 
atoms, reducing the accuracy of the detection of the applied force that induces the relative 
phase (p. In order to analyze the effect of statistical atom number fluctuations in the collisions 
of the solitons with the barrier, we simulate the soliton dynamics by incorporating the atom 
statistics in TWA. In Section 15.21 we show how the nonlinear soliton-barrier interaction 
increases the uncertainty in the relative atom number between the solitons and effectively 
cancels any potential benefit that could be obtained from the sensitivity of SN/N to variations 
of the phase (/>. 

4. Truncated Wigner Approximation 

In this section we outline the methods we employ for extending our study to include quantum 
fluctuations. A common approach in stochastic phase- space methods is to unravel the 
evolution dynamics into stochastic trajectories, each of which obeys the classical mean-field 
dynamics. Each trajectory is a representative of a probabilistic initial state distribution that 
is numerically generated by stochastic sampling. The probability distribution is selected 
so that one can synthesize as closely as necessary the quantum statistical correlations or 
thermal distribution of the atoms in the initial thermal equilibrium state. The phase-space 
representation that most accurately reproduces classical mean-field dynamics is the Wigner 
representation because of the 'correct amount' of quantum noise in the initial state p7l[36|. 

Here we introduce stochastic phase-space methods for the atom dynamics within TWA 
Il34l [35l [36l [37l [381 [39l l40l in order to investigate atom number statistics of the solitons, 
collapse of relative phase-coherence, and the condensate fragmentation. In our approach we 
follow the formalism of Ref. ||37|| . We consider a unitary evolution of the ID atomic gas 
without additional couplings to environment and replace quantum field operators (ij/, if/'^) by 
the classical fields {i//w, lAw) that are the stochastic phase-space representation of the full field 
operator. The nonlinear dynamics of each stochastic trajectory, sampled from the distribution 
of the initial conditions (we sample between 100 and 1000 stochastic realizations of i//w for 
each parameter regime), is then governed by the GPE evolution [Eq. ([I])]- The quantum 
fluctuations are stochastically included in the initial state ij/wix, t = 0). We solve the initial 
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state within the Bogoliubov approximation when the full quantum field operator is given by 

if/{x, 0) = 4fodo + ^ [ujix)aj - v*ix)a'[j . (25) 
j 

Here the phonon mode functions uj and Vj are calculated numerically. The stochastic initial 
state is obtained by replacing the phonon mode operators o-y, by classical stochastic 
variables fl37|] . Within TWA, ensemble averages and variances of observables correspond 
to quantum mechanical expectation values and uncertainties. Individual realizations of i^fw 
represent outcomes of individual experimental runs. The TWA method becomes more 
accurate when the modes are highly populated. However, the TWA simulations have been 
successful in providing qualitative descriptions of experiments [|48l on dissipative quantum 
dynamics in ID lattices even for significant ground state depletion and at very small atom 
numbers, as shown in Ref. 

Although the ID dynamic model neglects density fluctuations in the radial direction, it 
can also provide a good qualitative description of several experimental systems in which the 
radial confinement is not strong. This was demonstrated in Ref. [50J where experimentally 
observed atom number squeezing and reduced on-site atom number fluctuations were 
compared with numerical ID TWA simulation results. It has been suggested, however, that in 
a condensate collapse mechanism and soliton train formation a full 3D treatment of radially 
weakly confined systems can lead to distinct features in dynamics [fTTll . 

In order to calculate quantum mechanical expectation values and correlation functions 
for the soliton observables, we follow a method similar to that in Refs. [[37l ISTI and project 
the stochastic field ij/wix, t) onto some particular set of mode functions, here represented by 
the soliton solutions Tj{x, t). The amplitude of the jth soliton reads 

ajit) = ^ dx y*j{x, t)ifjw{x, t), (26) 

which corresponds to the stochastic Wigner representation for the annihilation operator of the 
soliton j. We define the normalized nonlinear mode functions T j{x, t) in such a way that the 
amplitudes a/0 incorporate information about the atom number and the phase of the soliton 
j 



g-i^j(r) \a\Ni r 1 

Yjix, t) = ilfjix, t) = i^^sech [|a|A^y(^ - q^)] (27) 



where ij/jix, t) is given by Eq. Q, and x, a, qj and Vj and are the variables defined in Sec. [31 
The phase of each soliton in each stochastic realization of the TWA dynamics is given by the 
complex phase of aj. In the classical limit, this phase corresponds to 0/t) in Eq. ([3]). 

In order to track the mode functions Tj{x, t) in time, we numerically calculate the soliton 
position in each realization of i/riy by integrating the centre-of-mass position J x\ifjw?'dx one 
or both sides of the barrier (if one or two solitons are present) and identify this with the 
soliton position(s), qj. We calculate the time derivative of qj along the each trajectory in ipw 
and identify this as the velocity Vj. The number of atoms Nj in each soliton is evaluated by 
integrating ji/'wP around the soliton's centre-of-mass position. 
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Quantum mechanical expectation values and correlation functions related to the soliton's 
atom number and the phase-coherence may then be calculated in the TWA simulations from 
the ensemble averages consisting of products of a/O's- The relative phase-coherence between 
two bright solitons j = 1,2 at times between the soliton splitting and recombination are 
obtained from the normalized first order correlation function 

c = \^ , (28) 

VW)(A^2) 

where the expectation value of the number of atoms in the jth soliton is obtained from 

<^,) = <a}a,)^-l/2. (29) 
The expectation value for the relative atom number 5N and its fluctuations are given by 

{6N) = (N2) -{Ni) = {ala2)w - {a\a,)w , (30) 

(A6Nf = {[N2 - Nif) - {N2 - Nif 

= {{a*2a2 - a*iaif)w - {a\a2 - a\a\% - l/l. (31) 

Within the Wigner representation, the stochastic variables correspond to symmetrically 
ordered expectation values, here denoted by (■ • that are obtained in the simulations as 
ensemble averages. The 1/2 term subtracted from the stochastic expectation value arises 
from the conversion between symmetrically and normally ordered operator averages. In our 
simulations typical occupation numbers for the solitons vary between 50 and 400000 atoms. 



5. Quantum dynamics 

We interpolate between regimes of significant quantum fluctuations and the classical mean- 
field limit, holding the nonlinearity \U\ [given by Eq. (fT6l) 1 constant and varying g/N (g can be 
adjusted in experiments, e.g., by varying the transverse confinement or the scattering length 
using a Feshbach resonance). In repulsively interacting atomic gases, g/N is proportional 
to the eff"ective interaction strength y = mg/ffn, where n denotes the ID atom density 
[f52|. that describes the strength of quantum fluctuations in the system. In tightly-confined 
attractively interacting atomic gases that we consider, \g/N\ also indicates the strength of 
quantum fluctuations. 

Fixing U to be constant fixes the number of excited-state atoms depleted from the ground 
state in the Bogoliubov approximation; hence, varying g/N is tantamount to varying the 
depleted atom fraction. When the magnitude of \U\ is large, i.e., the interactions between 
atoms are strong and the total atom number is small, quantum fluctuations are enhanced. By 
increasing \U\ we find significant eff'ects on the number fluctuations, phase decoherence and 
increased soliton position uncertainties for large \g/N\. The classical limit is N ^ 00, g ^ 0, 
when the depleted atom fraction becomes negligible. Here the dynamics is generally expected 
to approach the classical GPE description, the system becomes phase coherent, and the atom 
number statistics Poissonian. In the case of strongly nonlinear soliton-barrier interaction, 
however, we show that the nonhnear dynamics is very sensitive to any small perturbations 
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and the relative atom number statistics of the split solitons considerably exceed the classical 
shot noise in the limit N ^ oo, g ^ 0. 

5.1. Condensate fragmentation, phase-coherence and relative atom number statistics 

In Sec. |2]we discussed the importance of the relative phase-coherence between the solitons 
in the condensate collapse experiments leading to formation of soliton trains [111 [31 [H and in 
the experimental realization of an atom-optical soliton beam-splitter and interferometer that 
require a high level of coherence f3U . Here we study the effects of quantum fluctuations 
on the soliton phase-coherence and show how the coherence has important implications for 
the many-body entanglement in the system, and to possible quantum-enhanced precision 
measurements. We calculate the relative phase-coherence between the two solitons that are 
spatially separated by the barrier. We determine how well the phase-coherence is maintained 
during the spatial separation before the solitons are recombined and under which conditions 
the phase-coherence is lost and the two solitons represent the fragmentation of the initially 
phase-coherent condensate. 

Quantum fluctuations have been experimentally observed in BEC systems that are 
divided into spatially separated regions, e.g., by optical lattices [[531 HH [551 EH [571 EH [59l 
[60l[50l. The reduced tunneling amplitude for the atoms between adjacent sites and repulsive 
inter-atomic interactions lead to sub-shot-noise fluctuations in relative atom numbers between 
diff"erent lattice sites, along with reduced relative phase-coherence. Here we find that the 
interaction of the soliton with the laser barrier results in the splitting of the soliton into a 
reflected and transmitted soliton that exhibit non-classical relative atom number correlations 
even in the case of large atom numbers. We also examine how the relative atom number 
fluctuations between the solitons, when they exit the barrier the second time, are related to the 
decoherence rate of the solitons and show that the loss of the phase-coherence can be detected 
by monitoring the atom number statistics of the exiting solitons. 

In order to quantify the relative atom number fluctuations between the two split solitons 
we calculate the relative atom number squeezing obtained from 

f,^(^Mf^-H^. (32) 
4<Wj>(Af,> 

In the shot-noise limit, the relative atom number fluctuations follow a binomial distribution, 
and ^jj = I. Interactions between the atoms may suppress the atom number fluctuations such 
that < 1, in which case the fluctuations are number-squeezed. (When (^2) = (Ni), this 
condition is equivalent to {A6N)^ < N.) Alternatively, the nonlinear character of Eq. ([Tj) 
may cause the system dynamics to be sensitive to small perturbations, such that fluctuations 
are enhanced above the shot-noise limit (^^ > 1). We consider how depends on the 
strength of the nonlinearity by varying U [Eq. ([T6l) 1 in systems with fixed total atom number 
= 1200 and with centre-of-mass kinetic energy per atom Te - llhoj [defined in Eq. ([T71) 1 
as the solitons hit the barrier. We study the time evolution between the first two collisions 
of the solitons with the barrier, i.e., after a single soliton is split into two and before the 
first recombination at the barrier. We study barriers of diff"erent widths; as in the classical 
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simulations we use both a narrow barrier described by Eq. (fT5]) with width cr ^ 0.02 IZ, which 
closely represents a sharp ideal barrier, and a wider barrier with width cr = 0.14/ (close to the 
experimental values of Ref. [,3 111 ). Figure He) shows that for the lowest value of \U\ that we 
consider, the fluctuations reach the shot-noise limit (^^ = 1), as expected for a non-interacting 
condensate. As the interaction strength is increased to the values Sfio) < \U\ < 5fia>, we find 
the relative atom number fluctuations to be squeezed, with =^ 0.89(+0.04) at |f/| = 5fia>l 
when soliton is split by the narrow barrier (cr = 0.02 IZ). For the wider barrier case (cr = 0.14/), 
squeezing is weaker, with ^ 0.98(+0.05) for the same nonlinearity. We quote the statistical 
sampling errors in which occur due to the finite number of stochastic trajectories in each 
ensemble. For stronger interactions, the squeezing disappears, and becomes much greater 
than one as the interaction between the soliton and the barrier in the splitting process becomes 
more nonlinear and sensitive to small perturbations. It is possible that other splitting methods 
less sensitive to perturbations may produce relative atom number-squeezed solitons in the 
limit of large \ U\. 

From the TWA simulations we calculate the dynamics of the relative phase-coherence 
and the decoherence rate between the two split solitons by evaluating c [Eq. (|28]) 1 after the 
soliton splitting and before the re-collision of the two solitons at the barrier. The value of c 
close to one indicates a high degree of the relative phase-coherence and the decreasing values 
of c describe the loss of coherence and the condensate fragmentation. We also evaluate the 
spin-squeezing between the two solitons, associated with the fictitious pseudo-spin J = N/2 
of the two-state system of the two solitons. If (^A'^) = (e.g., between the soliton splitting and 
re-collision in the systems we consider), the spin-squeezing parameter is given by 

= (ASN)' , ■ ^ . (33) 

where 7^. = (alai + a\a2)ll and Jy = i{a\ai - a\a2)l2. The system is spin-squeezed with the 
relative atom number representing the squeezed quadrature if < 1, indicating a possibility 
for quantum-enhanced metrology [1611 |62l [631 [641 [65 J . The condition < 1 can also be 
identified as a signature of quantum many-particle entanglement between the two solitons 
[[66l. When both solitons have large populations, the phase-coherence is reflected in a large 
value of {JxY + (-A)^' such that small number fluctuations and high phase-coherence are 
tantamount to spin-squeezing [[57l . Since the relative phase between the reflected and the 
transmitted solitons 5(f) - n/l, the mean value of the pseudo-spin is aligned close to the y 
axis and - ^^/(sin50)^, where the phase-coherence (sin 50)^ could be experimentally 
measured by interfering the solitons or, as explained below, by monitoring the atom statistics 
of the output of the soliton beam-splitter after the second collision with the barrier. 

We calculate the relative phase-coherence c between the two solitons both in the case of 
a fixed nonlinearity \U\, when the initial quantum depletion is varied by changing the value 
of g/N, and in the case of a fixed total atom number A'', when we vary the nonlinearity \U\. 
We study the phase-coherence after splitting with the same barrier parameters with which we 
studied the relative atom number squeezing. The effect of the strength of quantum fluctuations 
on the decoherence rate of the two solitons is illustrated, for the case of the wide barrier, in 
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Figure 3. (a) Evolution of the relative phase-coherence c between two solitons, between the 
times (jL>t ^ 2 (immediately after splitting of the initial soliton by the Gaussian barrier of the 
same parameters as in Fig. [T] (e)-(i)) and <yf ^ 4 (immediately before the solitons re-collide 
at the barrier). The kinetic energy per atom as the solitons hit the barrier is the same as Fig. 
[Hand the nonlinearity \U\ = IShco. The different curves correspond to different values of the 
fluctuation parameter, corresponding to the atom numbers - 300 (magenta), - 600 (red), 
1200 (black), 8000 (blue) and 80000 (green). The fluctuation parameter g/N varies between 
approximately -1.7 x 10 '*/j<y/ and -2.3 x IQ^'^Hojl. (b) The phase-coherence between solitons, 
c [given by Eq. (l28Tll. as a function of nonlinearity \U\ (in the units of Hll)) with = 1200 at 
time ojt - 4, immediately before the solitons re-collide at the barrier The red line corresponds 
to a system split by the same barrier as in (a), and the black line to a system split by the 
narrower barrier used in Fig. [na)-(d). The barrier height is different for different values of 
\U\ chosen so as to split the soliton in two equal parts, (c) Relative atom number squeezing 
parameter [given by Eq. (l32T il measured at a time between the splitting of a soliton and re- 
collision for the same system as (b). The red and black lines are for the same wide and narrow 
barriers as in (b). The error bars show the statistical uncertainty in the value of arising due to 
the finite number of stochastic trajectories in each ensemble. The solitons in the narrow barrier 
case display stronger relative atom number squeezing and also exhibit spin-squeezing or non- 
classical entanglement, (d) The quantum uncertainty in the relative atom number between the 
solitons A6N, normalized by total atom number following the solitons' re-collision at the 
barrier as a function of their relative phase-coherence c at a time immediately before the re- 
collision, after phase imprinting by a potential given by Eq. jTli with amplitude V,,, ^ 3.4fiaj. 
The different points are for the same total atom number A^ but with the nonlinearity \U\ varied 
between the non-interacting and the strongly nonlinear systems. 
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Fig. [2]; a) for I i7| = 1 5 (the narrow barrier gives similar results). We show the time evolution 
of c between the initial splitting and the recombination of the solitons at the barrier, i.e., at 
the times 2 < ojt < 4. The different curves represent different depleted fractions of the atoms 
from the initial ground state of the single soliton, obtained by varying g/N. We find that 
high relative phase-coherence between the two split solitons is maintained during the entire 
evolution between the splitting and the recombination at the barrier for = 80000 (for this 
nonlinearity, the fluctuation parameter g/N - -2.3 x lO'^^hojl), even though the solitons are 
separated ~ 10/ apart. Here = 80000 represents almost the classical limit with c - 1 
and with a very slow decoherence rate. The other curves demonstrate a stronger decoherence 
rate and the condensate fragmentation, representing a transition from the classical to quantum 
limit as quantum fluctuations become progressively more prominent. In the case of strong 
quantum fluctuations {N = 600, such that g/N - -4.2 x \0~^fia)l) the coherence is already 
notably reduced immediately after the splitting (tot - 2) and decays rapidly. Moreover, we 
find that the decoherence rate is sensitive to the nonlinearity; for nonlinearities \U\ < 5fito 
the solitons remain phase-coherent with c > 0.9 even for large g/N (atom numbers as low 
as 600). We fit the curves in Fig. [Sfa) (where the wide barrier was used to split the solitons) 
to the function c(0 = A exp[-Yiojty] and obtain F ^ 1.6, 1.2, 0.22, and 0.03; a - 1.1, 1.6, 
1.6 and 1.3; A ^ 0.85, 0.89, 0.97 and 0.99 for the atom numbers = 300, 1200, 8000, and 
80000, respectively. For curves obtained using the narrow barrier, the corresponding values 
are F ^ 3.1, 1.5, 0.27, and 0.03; or ^ 2.1, 1.3, 1.4 and 1.3; A ^ 0.78, 0.89, 0.97 and 0.99. Here 
the fluctuation parameter g/N takes values between 1.7 x lO^'^hojl and -2.3 x 10~^ fiojl. 

We show in Fig. ^h) the relative phase-coherence c for a fixed total atom number 
A'^ = 1200, and kinetic energy per atom as the soliton collides with the barrier Te - llfioj, 
when we vary the nonlinearity \U\. The coherence is evaluated dX ojt ^ A - immediately 
before the solitons re-collide at the barrier - providing information about the coherence of 
the atoms at the soliton beam-splitter. Strong phase-coherence is retained for nonlinearities 
as strong as ^ Shco. When \ U\ is increased, the phase-coherence decreases monotonically 
and with a steep slope for \U\ > Sfioj (e.g., for the atom number A'^ = 1200, c < 0.2 for 
\U\ = IShoj). We find that the behaviour of the coherence c is quantitatively similar for both 
barrier widths. We fit the curve in Fig. Ob) corresponding to the narrow barrier to the function 
c(0 = Aexp[-F(|t/|/;ia;)"'] and obtain F ^ 1.5xlO-^a - 3.5 and A ^ 0.98. At\U\- SHcothe 
relative atom number fluctuations are most squeezed with ^ 0.89(+0.04) when split by the 
narrow barrier, and - 0.98(+0.05) when split by the wide barrier. Due to the high relative 
phase-coherence between the solitons we find that they also exhibit non-classical quantum 
correlations, or many-body quantum entanglement, as the system is also spin-squeezed when 
the narrow barrier is used, with ^ 0.95(+0.04). However, for the wider barrier case the 
spin-squeezing is lost, with ^ 1.08(+0.05). The example demonstrates how spin-squeezed 
solitons can be generated even with large atom numbers, providing a very different approach 
from the optical lattice experiments [571 EQI for producing spin-squeezed atomic samples. 

The loss of the relative phase-coherence between the two solitons could be 
experimentally measured by interfering the solitons. Here we propose an alternative method 
of experimentally measuring the loss of the relative phase-coherence that is directly based on 
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detecting the fluctuations in the outputs of the soliton beam-splitter. In the classical mean-field 
calculations we showed in Sec. [3] that the relative number of atoms in the solitons after the 
recombination at the barrier is dependent upon their relative phase, with the non-interacting 
limit given by Eq. (I24l) . By including quantum fluctuations in the stochastic TWA simulations 
we demonstrate here that the loss of phase-coherence is correlated with increased relative 
atom number fluctuations between the two solitons after the second collision with the barrier. 
Detecting the fluctuations in the relative soliton atom number could therefore provide a direct 
signature of the relative phase-coherence and the decoherence rate between the solitons. 

The large relative atom number fluctuations between the two solitons exiting the barrier 
are illustrated by the particular case shown in Fig. SI The relative phase-coherence between 
the solitons before they recombine at the barrier is in this case low with c < 0.2. We display in 
Fig.|4ta)-(b), and (c)-(d) two diff"erent examples of two individual stochastic realizations of i^fw 
with the same system parameters with the most of the atoms propagating in one realization 
to the negative x direction and in the other one to the positive x direction, indicating very 
large atom number fluctuations between different stochastic trajectories. The dependence 
of the relative atom number fluctuations after the second collision with the barrier on the 
phase decoherence is characterized more quantitatively by calculating quantum fluctuations 
from the ensemble averages of the stochastic simulations. We imprint the relative phase of 
6(f> - n/2 between the solitons as in Sec. [31 such that in the classical limit approximately two 
equal-sized solitons emerge from the barrier [and not only a single recombined soliton; see 
the non-interacting limit given by Eq. (|24l) 1. As shown in Fig. Od), the loss of the relative 
phase-coherence between the solitons is clearly signaled in an increase in the relative atom 
number fluctuations A6N/N after the second collision with the barrier. This provides a method 
of detecting the condensate fragmentation during the soliton splitting by the barrier. 

We show in Fig.[2a)-(b) the quantum mechanical expectation value of the relative atom 
number {6N) /2N in the solitons exiting the barrier (after the second collision of the atoms with 
the barrier in the recombination of the solitons) for different nonlinearities gN and total atom 
number A^, with the error bars displaying the corresponding quantum uncertainty {A6N)/2N. 
We also show in Fig. Oc) demonstrating how the relative atom number fluctuations 
typically exceed the shot-noise limit by an order of magnitude. As well as coinciding with 
increased fluctuations in the relative atom number of the outgoing solitons, the loss of phase- 
coherence also coincides with a shift in the expectation values {6N)/2 towards zero, since 
due to the negative curvature ^A'^ as a function of 6(f>, a symmetric distribution of phases 6^ in 
different stochastic realizations within TWA always lead to a smaller value of {6N). This is 
illustrated by the points with large error bars in Fig. [5ta)-(b) which lie along a flatter profile 
than the curves for smaller quantum fluctuations. 

5.2. Quantum Bright Soliton Interferometer 

In Sec. |3]we described the classical behavior of a soliton interferometer where during their 
spatial separation, the solitons can experience a relative phase produced by a small force, 
e.g., gravity. Such a scheme can be implemented by aligning the axial direction of our trap 
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Figure 4. Individual stochastic realizations of the Wigner density liAiyp within the TWA 
representing possible outcomes of individual experimental runs of the soliton-barrier collision 
dynamics for mean atom number = 600 and nonlinearity \U\ - IStiu, with (a)-(b) a narrow 
barrier of the same dimensions as in Fig. [na)-(d), and in (c)-(d) a wider barrier of the same 
dimensions as in Fig. [ne)-(i)- Although the physical parameters in (a)-(b) and also in (c)-(d) 
are equal, the trajectories correspond to the situations where the solitons exit in the opposite 
directions from the barrier due to quantum fluctuations. 

orthogonal to the direction of the earth's gravitational field, such that the solitons are separated 
by a vertical displacement 6x between times and tf. The solitons experience a potential 
difference of amplitude V,,, = mg6x which classically imprints the phase (p ~ mg6x(tf - ti)l%, 
where g is the local acceleration due to gravity. Outside the classical regime, this phase will 
be subject to quantum fluctuations resulting in the phase decoherence described in Sec. 15. 1[ 
providing one limit on the interferometric measurement. The accuracy will be limited also 
by the relative atom number fluctuations, described by the relative atom number squeezing 
parameter in Eq. (|32l) which occur due to the sensitivity of the soliton-barrier collision to 
small fluctuations. In Sec. [3] we showed steepening of the gradient of the curve 5N IN with 
increasing nonlinearities \U\ for applied phases (p below a maximum value. Here we show 
any benefit to the sensitivity of the phase measurement due to this steep gradient is eff"ectively 
cancelled by increased number fluctuations in the outgoing solitons. 
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Figure 5. (a)-(b) Relative atom number between outgoing solitons 6N normalized by twice 
the total atom number 2N after collision at the barrier in the soliton recombination following a 
phase-imprinting by a potential given by Eq. (|2TI) with maximum value y„, (in the units of hco). 
In (a) the barrier is narrow, with the same width as in Fig[na)-(d), and in (b) the barrier has the 
same width as in Fig.[TJe)-(i). The lines are from GPE simulations: the nonlinearity takes the 
values: |f/| = (black line), 3/i<w (red line) and lOhuj (blue line). The points and error bars of 
corresponding colours indicate the quantum mechanical expectation values and uncertainties 
A6N/2N for the TWA simulations of the same nonlinearity and total atom number = 1200. 
Note that the data for \ U\ = lOhco do not sit on a curve. Additional points with cyan error bars 
are for nonlinearity \U\ = lOfiw and total atom number = 800000, which sit on the curves 
produced in the GPE. (c) Relative atom number squeezing parameter as a function of total 
atom number for outgoing solitons. The applied phases are imprinted by the potential Eq. 
(ISTT l with y„, = (full line), Vm - 3.4fKD (dashed line), bStiu (dot-dashed line) and lO^ai 
(dotted line). The barrier width is the same as in (b), and the nonlinearity \U\ - 3fia), such that 
the system has the parameters of the red curve in (b). 

In the limit of negligible phase fluctuations, an interferometer's accuracy, A0, is given 
in terms of the small phase differences it can resolve. The standard quantum limit for a 
phase measurement by a classical interferometer with particles is = l/yfN [i65l . 
One may apply a controlled potential to one soliton, such that the relative phase between 
solitons has some particular expectation value (50) and evaluate the ability of the scheme 
to distinguish two relative phases: (6^) and {6^i). We consider regimes of high phase- 
coherence, where the uncertainty is dominated by the enhancement of fluctuations during 
the collision with the barrier: even for very coherent solitons, is significant. Thus, we 
may characterize the accuracy of the interferometer by mapping 6N as a function of (6^) 
and defining the closest resolvable atom numbers (^A'^) and {6Ni) whose separation must 
be of the order of the fluctuations \{6N) - {6Ni)\ - A6N. We then define the accuracy 
as = ((501 ) - (50), where (50), (50i) are the relative phases that would produce 
the expectations in the relative atom numbers {6N), {6Ni) in the output solitons. We 
approximately obtain A0 ^ A6N/ ^d{6N)/d{6^)^ in which case the accuracy can be improved 
by having a very sensitive dependence of the relative atom number on the induced relative 
phase or by suppressed relative atom number fluctuations. 

In the non-interacting (U = 0) case, there is no phase diffusion during the soliton 
dynamics and no nonlinear interaction between the solitons and the barrier, and the output 
will be limited by shot noise such that the relative atom number fluctuates according to the 
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classical binomial statistics with = 1 . For this system we consider a force which in classical 
GPE simulations would imprint a phase close to n/2, and within TWA induces the quantum 
expectation of the relative phase {6(f)) = (which contains canceling contributions from the 
phase shift of the reflected soli ton and from the phase imprinted by the potential). Here the 
gradient of 6N as a function of (Sep) is largest, which should give the best accuracy for the 
resolvable phase diff"erence A0. One can evaluate Acp by mapping the fluctuations in the atom 
number and the gradient of the curve 6N((p) [e.g., in Fig. [3a)-(b)] to the phase uncertainty 
by the formula A0 ^ A6N/ [d{6N)/d6(()]. In the non-interacting case this takes the form 
A0 ^ 1 / VA^ (the standard quantum limit) in the limit where is large. 

In Sec. [3] we showed classical GPE simulation results for interacting solitons that 
demonstrated high sensitivity of 6N on 6(() with large values of d{6N)/d(6(f>). The classical 
GPE does not provide any information about the atom statistics and if we assume a classical 
coherent state with A6N - VA^, we would expect an interferometer accuracy that is better 
than in the ideal gas case or the standard quantum limit. Our stochastic TWA simulations, 
however, show that A^A'^ is substantially increased due to nonlinear interaction of the solitons 
with the barrier whenever d{6N)/d{6^) is large. This happens even in the classical limit of 
very large atom numbers and negligible phase decoherence with c - 1 . 

In the TWA simulations the phase decoherence is not always negligible and we study 
the regimes that we found in Sec. 15. II to have strong phase-coherence between the solitons 
throughout their evolution. This is the case either for weak interactions {\U\ < 5fico) for all 
atom numbers we considered A'^ > 600, or in the classical limit of large A^ for any gN. Within 
the regimes of high phase-coherence, we consider situations where d{6N)/d{6^) is large, 
corresponding to large classical values of d(6N) / d(6(f>) found in Sec. [31 (Here the interactions 
cause the shape of [6N/2N] to distort, such that its maximum occurs at smaller applied 
forces, apart from in the exceptional case of a large nonlinearity and wide barrier, where 
the curve is also shifted with respect to the applied force.) If were fluctuating according to 
the classical binomial distribution, the increase in d{5N) I d{54>) would decrease the phase 
uncertainty below the standard quantum limit. However, the interactions also afi'ect the 
number fluctuations [Fig.[2a)-(b)] and we find from TWA simulations that when the gradient 
is steepest, the relative atom number fluctuations are largest, in fact, within these regimes we 
do not find suppressed fluctuations in the output solitons for any applied force; sensitivity 
of the re-collision of the solitons at the barrier to small perturbations leads to increased 
fluctuations in the outgoing solitons for any atom number. For instance, for \U\ = 3fioj, in 
Fig- He) we find that for the finite imprinted phases, even in the classical limit of large atom 
number A^^, for fixed U, converges to approximately 36 for the wider barrier (^^ converges 
to 40 for a narrow barrier) at the applied potential for which the number fluctuations have the 
largest gradient (here {5N)/2N - 0). From A(f> ^ A6N/ [d{SN)/d{6^)], we obtain an optimum 
phase accuracy A(f> of approximately 5 times the standard quantum limit. Interestingly, when 
no potential is applied (V,„ = 0) and we have {6N) - -N, for a narrow barrier the fluctuations 
reflect the initial fluctuations in the total atom number in the classical limit N ^ oo and U held 
constant. For a wider barrier, the fluctuations remain finite but smaller than for those in the 
systems with finite V,„. However, in this limit, there is only negligible relative population in the 
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outgoing left soliton, and the relative atom number fluctuations are hard to track numerically. 
We consider the more strongly nonlinear system with \U\ = lOfiaj and A'^ = 800000; here 
there is approximately a three-fold increase in the gradient of (SN) as a function of (0) from 
the non-interacting case. However, strong number fluctuations (e.g., - 110 for the narrow 
barrier) destroy the ability to measure phase fluctuations below the standard quantum limit 
despite the increase in the gradient of {6N). Where the gradient of Fig. [5ta) is maximal, 
this would result in a phase uncertainty nearly 7 times the standard quantum limit (which 
would occur if the relative atom number fluctuations were binomial and the gradient was not 
steepened with respect to the non-interacting case). This phase uncertainty is also 21 times the 
value which would occur if the system had binomial number fluctuations accompanied by the 
steep gradient. In simulations extending over multiple periods of the soliton dynamics in the 
harmonic trap, we find that the number fluctuations increase with repeated soliton collisions 
with the barrier, which we attribute to the nonlinear interactions between the solitons and the 
barrier. 

5.3. Sub-shot-noise soliton interferometer 

Sub- shot-noise atom number fluctuations provide a central tool for quantum interferometric 
sensors that can operate in precision measurements at higher accuracies than the standard 
quantum limit of classical atom interferometers ll62l [63l [64l l65l . In particular, the phase 
measurement errors may reach the Heisenberg limit, for the error scaling as 1/A^ with the 
occupation number A'^, as compared to the standard shot-noise error l/N^^^. Reduced atom 
number fluctuations find interferometric applications in atomic systems ll6Tl[67l and the spin- 
squeezing of a soliton pair (Sec. [5]) suggests a possibility for an atomic soliton interferometer 
realization resulting in quantum statistical properties that can overcome the shot-noise limit 
of classical atomic samples. Although the soliton interferometer described in the previous 
section cannot achieve sub-shot-noise accuracies, we propose a scheme that potentially could 
provide a soliton interferometer overcoming the standard quantum limit that is based on 
analogous state manipulations as Ref. |f6T|. After producing a spin-squeezed soliton pair 
(Sec. [S]), one may perform interferometric manipulation of the internal states of the atoms, 
instead of letting the solitons collide at the barrier. Using a n pulse one may transfer the 
atoms in one of the solitons to a difl^erent internal state; then, by removing the laser barrier, 
let the solitons overlap as they reach the centre of the trap. During this spatial overlap 
one may turn up a Raman coupling between the two internal states for time t, inducing a 
rotation in the Bloch sphere by exp(z'rA:7;i) evolution. In order to minimize the efl^ects of 
interspecies interactions between the two solitons during the overlap, we require the strength 
of the Raman coupling k to be sufficiently strong compared to the interaction strength between 
the two states. By choosing tk = n/2, one can thus transfer number squeezed states to phase 
squeezed states, such that any accumulated relative phase can be measured with accuracy 
below classical precision after another n/2 rotation of the Bloch sphere. 
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Figure 6. (a)-(c) Quantum uncertainty in soliton momentum Ap as a function of total atom 
number during evolution in a harmonic potential (with no Gaussian barrier included) at 
times ojt ^ (a) 0.15, (b) 0.70 and (c) 1.3. (d) Quantum uncertainty in soliton position Ax at 
time tot ^ 0.15 (Ax remains approximately the same at later times). The soliton has kinetic 
energy per atom as it reaches the trap centre Te - llfio) and nonhnearity \ U\ - 3hcj. Ap and 
Ax reach their classical (GPE) limits, represented by the dotted lines, before = 80000 - well 
outside the plotting range. We found the statistical errors on these data to be negligibly small. 
For example, for A' = 600, at ^ 1 .3 the error in {Ap)l/ti is approximately 1 .7 x 10"^ and the 
error in Ax/l is approximately 2.1 x 10"^. 



5.4. Quantum statistics of the soliton's position and momentum 

In the previous sections we studied quantum effects on the relative atom number fluctuations 
between the solitons and on the relative phase-coherence in the soliton-barrier system. Here 
we investigate the effects of quantum fluctuations on soliton position and momentum by 
extracting the quantum expectation values and uncertainties for each soliton's position and 
momentum from TWA. As described in Sec. |4] we numerically track the centre-of-mass 
position qj and velocity Vj of each soliton, as well as the atom number Nj in every soliton, as a 
function of time in every stochastic realization. The quantum mechanical expectation values 
of the first and second powers of soliton position and momentum: {xj), (xj), (pj) and (pj) may 
then be obtained by calculating ensemble averages of (qj), (q^j), (vj) = (dqj/dt), etc., from 
the TWA simulations and using the relationships Eqs. ([8l)- (fT2)) over many realizations, so that 



Quantum dynamics of atomic bright solitons 



24 




Figure 7. Trajectories of the quantum mechanical expectation values and uncertainties in 
soliton position (a) and momentum (b) for a soliton of the same kinetic energy and split by a 
barrier of the same parameters as in Fig.[T]into two solitons which subsequently re-combine 
at the barrier. The quantum expectation value of the position in (a) and the momentum in 
(b) are marked by the dark solid lines, whereas the quantum uncertainties are defined by 
the corresponding shaded areas. The nonlinearity \U\ - TiTkjj and the total number of atoms 
= 1200 (magenta) and = 300 (yellow) such that ^/A? = -2.8xlO"^Sw/ and -S.SxlO^'^^iw/ 
respectively. 



the quantum mechanical position and momentum uncertainties for the solitons approximately 
satisfy 

{Apjf ^ {p]) 

where A^y = ((^/) - iqjY) is the uncertainty of the centre-of-mass position of the soliton 
obtained from the stochastic TWA simulations. The actual values of the uncertainties may 
slightly differ from the approximate right-hand-side expressions of Eqs. (l34l) and (l35l) due 
to the deformations of the shape of the soliton wavefunction as a consequence of quantum 
fluctuations and nonlinear interactions. The position and momentum uncertainties therefore 
contain contributions from the width of the soliton's wavefunction as well as the quantum 
fluctuations in each soliton's centre-of-mass position and velocity. In the GPE limit, {lS.Xjf = 
(Axjf, = ijn^/iUa^Nj) and (Apjf = {Apj)] = nVNj/i3l% as given by Eqs. ^ and ([H). 
Hence, the product ApjAxj is larger than the Heisenberg uncertainty limit by a factor of n/3. 
We find that in the TWA simulations quantum fluctuations increase this product further above 
the Heisenberg uncertainty limit, e.g., for \U\ = Shco and for the small atom number = 50, 
the product ApjAxj is approximately 6% larger than in the GPE limit. 

We first consider the periodic dynamics of a single soliton. We investigate the effect 
of quantum fluctuations by varying g/N for a fixed U, so that the initial states of the TWA 



Ifa'^N^ 

(Pjf^^^+miAvjf, (35) 



Quantum dynamics of atomic bright solitons 25 




-5 5 -5 5 




{p)i/h,iAp)i/n {p)i/h,{Ap)i/h 

Figure 8. (a)-(b) Individual stochastic realizations of the Wigner density li/'ivP within the 
TWA for nonlinearity \U\ = IShoj, in the regimes where solitons travel in different directions 
after colliding with a barrier described by Eq. ( fTSl ). The barrier height Vq ^ 6.Qftai and width 
cr - 0.32/, the total atom number = 600 and the kinetic energy per atom of the soliton as 
it hits the barrier Te - S.lhcj. (c)-(f) Quantum expectation values and uncertainties in soliton 
position (c)-(d) and momentum (e)-(f) for the same nonlinearity U and potential barrier as 
(a)-(b), with total atom number values = 600 in (c),(e) and 8000 in (d),(f). The quantum 
expectation value of the position in (c),(d) and the momentum in (e),(f) are marked by the dark 
solid lines, whereas the quantum uncertainties are defined by the corresponding shaded areas. 
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Figure 9. Trajectories for the same parameters as Fig. [8] but with atom numbers = 1200 in 
(a),(c) and = 2000 in (b),(d). The quantum expectation value of the position in (a),(b) and 
the momentum in (c),(d) are marked by the dark solid lines, whereas the quantum uncertainties 
are defined by the corresponding shaded areas. 

simulations range from the classical limit of a negligible fraction of depleted atom number 
to the cases in which quantum fluctuations have a notable influence on the position and 
momentum uncertainties. In the quantum limit of large |^/A^|, the quantum uncertainty in 
soliton position Ax is large, and decreases with increasing A^, reaching the GPE limit before 
= 80000 for the studied values of U. Surprisingly, for very large \g/N\, at early times 
in the soliton dynamics, the uncertainty in velocity is smaller than in the classical mean- 
field limit [illustrated in Fig. |6ta) for \U\ = 3fito with A'^ varying between 100 and 2000]. 
The initial velocity uncertainty in this case is dominated by the width of the classical soliton 
wavefunction solution in momentum space [i.e., the first term in Eq. (|35l) 1. and at small A'^ 
there are relatively fewer atoms in each soliton due to quantum fluctuations, thus decreasing 
the nonlinear effect of interactions to broaden the velocity distribution. During the soliton 
evolution in the trap, the momentum uncertainty is aff"ected by the dissipation of soliton 
velocities within the ensemble, and Apj become large in the quantum limit of small A'^ 
(when \ U\ is constant) [see Fig.[6ta)-(c)]. A similar eff"ect was found for trapped dark atomic 
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solitons 10911681 . 

Quantum fluctuations also appear to have an unexpected effect on the quantum 
expectation values of the soliton momentum [Fig.|7]^b)], with (pj) decreasing in the quantum 
limit of smaU A^. For instance, at = SHco and = 300 (g/N = -3.3 x IQ-^fitol), the 
maximum expectation value of the soliton's speed is reduced by a few percent compared 
with the = 1200 (g/N = -2.08 x W'^hcol) case of the same nonlinearity \U\ = 3fia). 
For dark solitons in tightly-confined bosonic atomic gases the speed is reduced due to phase 
fluctuations across the phase kink and a nonlinear relationship between the value of the phase 
and the speed [|39l[68l . In the case of a bright soliton the changes in the quantum expectation 
value of the momentum, however, are more likely to be an indication of dissipative centre-of- 
mass transport of the atoms [i48ll49l . 

Previous quantum mechanical treatments of bright solitons with very small atom 
numbers interacting with a barrier have found Schrodinger-cat states of solitons, i.e., quantum 
superpositions of solitons traveling in both directions after hitting the barrier [|25ll26l . Such 
superposition states cannot be represented by classical stochastic phase-space methods. In 
systems with large atom numbers which we consider in this work, however, any interaction 
with the environment and a loss of atoms from the system would rapidly decohere a 
macroscopic superposition state into a statistical mixture of diff"erent atom number states 
[|69l iTOl that can be represented by TWA. In the simulations we find regimes in which the 
soliton propagates either left or right depending on the particular stochastic realization of i//w 
even though all the physical parameters remain the same [Fig.[8ta)-(b)]. We consider a barrier 
given by Eq. (fTSi) with height Vo - 6.0fito and width cr = 0.32Z, wider than the barriers 
previously considered in this paper, and approximately double the width of the barriers in the 
recent experiment [|3TI . We use soliton displacements of Xd - 3.21 such that the kinetic energy 
of the soliton as it hits the barrier is Te - S.lhoj. For small values of \U\, the barrier splits 
a soliton into two in each realization of i/^vk- However, increasing the nonlinearity such that 
\U\ = IShoj, the solitons are either entirely transmitted, or entirely reflected from run to run in 
different realizations of tffw This is indicated in the increased quantum uncertainty in soliton 
position after the soliton hits the barrier [Fig. Oc)]. By increasing N, while at the same 
time keeping U constant, we find the classical limit where the soliton is reflected in every 
realization of tffw and the position and momentum uncertainties remain constant during the 
soliton dynamics; we find that this limit is reached at = 80000. The momentum uncertainty 
Api is dominated by the width of the classical solution to the soliton wavefunction, rather than 
by the dissipation of the soliton momenta within the ensemble. However, fluctuations have 
a marked effect on the momentum expectation values, deforming the oscillatory behavior of 
(p) [see Figs SeXf)]. 

6. Concluding remarks 

We studied a physical realization of an atom-optical soliton beam-splitter and an 
interferometer in a harmonic trap where a repulsive laser barrier splits a soliton into two 
and recombines the split solitons. We found that the classical mean-field analysis based on 
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the GPE simulations, that does not provide any information about the atom statistics of the 
solitons, is not sufficient to estimate the accuracy of the interferometer even in the limit of a 
very weak quantum depletion and phase decoherence. This is because the nonlinear soliton- 
barrier interactions lead to super-Poissonian atom number fluctuations that substantially 
deteriorate the performance of the interferometer. The result may have implications also to 
other possible interferometer applications in nonlinear systems in which a high sensitivity 
of the interferometer output on the induced phase difl^erence suggests a possibility for high 
precision measurements, emphasizing the importance of a detailed statistical analysis of the 
particle number fluctuations beyond the mean-field theory. 

We also studied the periodic dynamics of oscillating bright soliton in the regime of 
considerable quantum fluctuations. Quantum fluctuations were found to have a notable efl'ect 
on the splitting, dynamics and recombination of the solitons, resulting, e.g., in the condensate 
fragmentation and in non-classical correlations (spin-squeezing) between the solitons even 
for large atom numbers. The bright soliton dynamics also exhibits large quantum fluctuations 
of the soliton's position and velocity. At the same time the soliton can represent a truly 
mesoscopic quantum object that may consist of hundreds or even thousands of atoms. This 
provides an interesting opportunity to vary the size of the solitons and use them as an 
experimental probe of the quantum-classical interface [TTTII . 

Acknowledgements 

We acknowledge financial support from the Leverhulme Trust and the EPSRC and discussions 
with Randall G. Hulet. 

References 

[1] Haus H A and Wong W S 1996 Rev. Mod. Phys. 68 423 
[2] Stegeman G I and Segev M 1999 Science 286 1518 

[3] Khaykovich L, Schreck F, Ferrari G, Bourdel T, Cubizolles J, Carr L D, Castin Y and Salomon C 2002 
Science 296 1290 

[4] Sti-ecker K E, Partridge G B, Truscott A G and Hulet R G 2002 Nature 417 150 

[5] Cornish S L, Thompson S T and Wieman C E 2006 Phys. Rev. Lett. 96 170401 

[6] AbduUaev F K, Gammal A, Kamchatnov A M, and Tomio L 2005 Int. J. Mod Phys B 19 3415 

[7] Eiermann B, Anker Th, Albiez M, Taglieber M, Treutlein P, Marzlin K-P and Oberthaler M K 2004 Phys. 

Rev. Lett. 92 230401 
[8] Dutton Z, Budde M, Slowe C and Hau L V 2001 Science 293 663 
[9] Busch Th and AngUn J R 2001 Phys. Rev Lett. 87 010401 
[10] Anderson B R Haljan P C, Regal C A, Feder D L, ColUns L A, Clark C W and Cornell E A 2001 Phys. 

Rev Lett. 86 2926 

[11] Becker C, Stellmer S, Soltan-Panahi R, Dorscher S, Baumert M, Richter E M, Kronjager J, Bongs K and 

Sengstock K 2008 Nature Phys. 4 496 
[12] Shrestha U, Javanainen J, and Ruostekoski J 2009 Phys. Rev Lett. 103 190401 

[13] Yan D, Chang J J, HamnerC, Kevrekidis P G, Engels P, Achilleos V, Frantzeskakis D J, Carretero-Gonzalez 

R and Schmelcher P larXivrt 10443591 
[14] Buljan H, Segev M, and Vardi A 2005 Phys. Rev. Lett. 95 180401 

[15] Al Khawaja U, Stoof H T C, Hulet R G, Strecker K E, and Partridge G B 2002 Phys. Rev Lett. 89 200404 



Quantum dynamics of atomic bright solitons 



29 



[16] Sti-eltsov A I, Alon O E and Cederbaum L S 201 1 Phys. Rev. Lett. 106 240401 

[17] Dabrowska-Wuster B J, Wuster S and Davis M J 2009 New J. Phys. 11 053017 

[18] Lai Y and Haus H A 1989 Phys. Rev. A 40 844 

[19] Lai Y and Haus H A 1989 Phys. Rev A 40 854 

[20] Dmmmond P D and Carter S J 1987 J. Opt. Soc. Am. B 4 1565 

[21] Werner M J and Drummond P D 1997 Phys. Rev A 56 1508 

[22] Fini J M, Hagelstein P L and Haus H A 1998 Phys. Rev A 57 4842 

[23] Corney J F and Drummond P D 2001 / Opt. Soc. Am. B 18 153 

[24] Jenkins S D, Salerno D, Minardi S, Tamosauskas G, Kennedy TAB and Di Trapani P 2005 Phys. Rev. 
Lett. 95 203902 

[25] Streltsov A I, Alon O E and Cederbaum L S 2009 Phys. Rev A 80 043616 

[26] Weiss C and Castin Y 2009 Phys. Rev Lett. 102 010403 

[27] Click J A and Can- L D, arXiv: 1 1 05 .5 1 64 

[28] Lewenstein M and Malomed B A 2009 New J. Phys. 11 1 13014 

[29] Ahufinger V and Sanpera A, 2005 Phys. Rev Lett. 94 130403 

[30] Lee R-K, Ostrovskaya E A, Kivshar Y S and Yinchieh L 2005 Phys. Rev A 72 033607 
[31] Dyke P, Sidong L, Pollack S, Hulet R G, unpublished 

[32] Cornish S L, Parker N G, Martin A M, Judd T E, Scott R G, Fromhold T M and Adams C S 2009 Physica 
D 238 1299 

[33] Parker N G, Martin A M, Cornish S L and Adams C S 2008 J. Phys. B 41 045303 
[34] Drummond P D and Hardman A D 1993 Europhys. Lett. 21 279 

[35] Steel M J, Olsen M K, Plimak L I, Drummond P D, Tan S M, Collett M J, Walls D F and Graham R 1998 

Phys. Rev A 58 4824 
[36] Sinatra A, Lobo C and Castin Y 2002 J. Phys. B 35 3599 
[37] Isella L and Ruostekoski J 2006 Phys. Rev A 74 063625 

[38] Blakie P B, Bradley A S, Davis M J, Ballagh R J and Gardiner C W 2008 Adv in Phys. 57 363 
[39] Martin A D and Ruostekoski J 2010 Phys. Rev Lett. 104 194102 
[40] Polkovnikov A 2010 Ann. Phys. (NY.) 325 1790 
[41] Al Khawaja U and Stoof H T C, 201 1 New J. Phys. 13 085003 
[42] Billam T P, Cornish S L and Gardiner S A 201 1 Phys. Rev A 83 041602(R) 
[43] Zakharov V E and Shabat A B 1971 Zh. Eksp. Tear. Fiz. 61 1 18 [1972 Sov Phys. JETP 34 62] 
[44] Faddeev L D and Takhtajan L A 1987 Hamiltonian Methods in the Theory of Solitons (Berlin Heidelberg: 
Springer- Verlag) 

[45] Martin A D, Adams C S and Gardiner S A 2008 Phys. Rev A 77 013620 
[46] V V Efimov, O V Ivanov and D I Sementsov 2001 J. Opt. A 3 514 
[47] Gardiner C W and Zoller P 1991 Quantum Noise (Berhn Heidelberg: Springer- Verlag) 
[48] Fertig C D, OHara K M, Huckans J H, Rolston S L, Phillips W D and Porto J V 2005 Phys. Rev Lett. 
94 120403 

[49] Ruostekoski J and Isella L 2005 Phys. Rev Lett. 95 1 10403 

[50] Gross C, Esteve J, Oberthaler M K, Martin A D and Ruostekoski J 201 1 Phys. Rev A 71 09040 1(R) 
[51] Isella L and Ruostekoski J 2005 Phys. Rev A 72 01 1601 (R) 

[52] Kheruntsyan K V, Gangardt D M, Drummond P D and Shlyapnikov G V 2003 Phys. Rev Lett. 91 040403 
[53] Orzel C, Tuchman A K, Fenselau M L, Yasuda M, and Kasevich M A 2001 Science 291 2386 
[54] Gerbier F Foiling S, Widera A, Mandel O, and Bloch I 2006 Phys. Rev Lett 96 053615 
[55] Li W Tuchman A K, Chien H-C and Kasevich M A 2007 Phys. Rev. Lett 98 040402 
[56] Sebby-Strabley J, Brown B L, AnderUni M , Lee P J, Phillips W D, Porto J V and Johnson P R 2007 Phys. 
Rev Lett 98 200405 

[57] Esteve J, Gross C, Weller A, Giovanazzi S and Oberthaler M K 2008 Nature 455 1216 
[58] Gemelke N, Zhang X, Hung C-L and Chin C 2009 Nature 460 995 

[59] Bakr W S, Peng A, Tai M E, Ma R, Simon J, Gillen J I, Foiling S, PoUet L and Greiner M 2009 Science 
329 5991 



Quantum dynamics of atomic bright solitons 

[601 Sherson J F , Weitenberg C, Endres M, Cheneau M, Bloch I and Kuhr S 2010 Nature 467 68 

[61] Gross C, Zibold T, Nicklas E, Esteve J and Oberthaler M K 2010 Nature 464 1 165 

[62] Wineland D J, Bollinger J J, Itano W M and Heinzen D J 1994 Phys. Rev. A 50 67 

[63] Holland M J and Burnett K 1993 Phys. Rev. Lett. 71 1355 

[64] Bouyer P and Kasevich M A 1997 Phys. Rev. A 56 R1083 

[65] Giovannetti V, Lloyd S and Maccone L 2004 Science 306 1330 

[66] S0rensen A, Duan L-M, Cirac J I and ZoUer P 2001 Nature 409 63 

[67] Riedel M F, Bohi, Li Y, Hansch T W, Sinatra A and Treudein P 2010 Nattire (London) 464 1 170 
[68] Martin A D and Ruostekoski J 2010 New J. Phys. 12 055018 
[69] Zurek W H, Phys. Today 1991 44 (10) 36, and references therein. 
[70] Walls D F and Milbum G J 1985 Phys. Rev A 31 2403 
[71] Anglin J 2008 Nature Phys. 4 437 



